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For neutrons bound inside nuclei, baryon instability can manifest itself as a decay into unde- 
tectable particles (e.g., n vvv ), i.e., as a disappearance of a neutron from its nuclear state. 
If electric charge is conserved, a similar disappearance is impossible for a proton. The existing 
experimental lifetime limit for neutron disappearance is 4-7 orders of magnitude lower than the 
lifetime limits with detectable nucleon decay products in the final state Q. In this paper we 
calculated the spectrum of nuclear de-excitations that would result from the disappearance of 
a neutron or two neutrons from ^^C. We found that some de-excitation modes have signatures 
that are advantageous for detection in the modern high-mass, low-background, and low-threshold 
underground detectors, where neutron disappearance would result in a characteristic sequence of 
time- and space-correlated events. Thus, in the KamLAND detector a time-correlated triple 
coincidence of a prompt signal, a captured neutron, and a /?''' decay of the residual nucleus, all 
originating from the same point in the detector, will be a unique signal of neutron disappearance 
allowing searches for baryon instability with sensitivity 3-4 orders of magnitude beyond the present 
experimental limits. 

PACS number(s): 24.80.-Fy, 13.30.Ce, 14.20.-c, 24.10.Pa 



I. INTRODUCTION 



Baryon asymmetry of the universe Q and the idea of 
unification of particle and forces 0, are the two global 
concepts that motivated the nucleon decay experimental 
searches 6] for more than two decades. The particular 
modes of nucleon decay are not known a priori, although 
some are preferred by theoretical models |j, la, H 01 • The 
results from several generations of nucleon decay search 
experiments 0, including the most spectacular recent 
Super-Kamiokande results @, have set stringent limits 
on many exclusive nucleon decay modes. For example, 
the lifetime limit for the decay mode p ^ e+ + 7r°, estab- 
lished by the Super-K Collaboration, reaches > 5 • 10'^^ 
years 10]. However, the mode-independent nucleon life- 
time limit |1| is established experimentally only at the 
level of > 1.6-10^^ years [ij, which is more than 8 orders 
of magnitude poorer than the previous limit. That indi- 
cates that experimental searches for certain decay modes 
need to be extended. 

The Particle Data Group identifies 75 possible modes 
of nucleon decay (baryon instability) H that respect the 
conservation laws of electric charge, energy-momentum, 
and angular momentum. The lifetime limits for most of 
these modes lie above lO'^" years, but there are a few 
significant exceptions. These are the modes correspond- 
ing to the decay of one neutron or two neutrons into 
neutrinos: n — » "iv, n — » bv, and nn 2v. Lifetime 
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limits for these modes are several orders of magnitude 
lower P|. Generally speaking, these modes can be inter- 
preted as disappearance of one or two neutrons from their 
intra-nucleus state into any undetectable particles. Fol- 
lowing the spirit of the paper of Evans and Steinberg [TJ] 
and using the assumptions of the conservation laws men- 
tioned above, one could think that the present "mode- 
independent" nucleon decay limit is determined by 
these neutron disappearance modes. We believe that ex- 
perimental searches for these modes and the improve- 
ment of mode-independent nucleon lifetime limit are as 
fundamentally important as searches for particular ex- 
clusive nucleon decay modes with super-giant detectors. 

Moreover, decays like N! — > lepton+lepton+antilepton 
arise m certain unification schemes Sill as a result of 
(B — L) non-conserving interactions that might be more 
natural for the explanation of baryon asymmetry of uni- 
verse than {B — L)-conserving transitions. The question 
of whether {B ^ L) is violated in nature and at what 
energy scale can be answered only by experiments |l3l |. 
Nucleon decay into five and more leptons was theoreti- 
cally considered in f3|; dinucleon decay into two leptons 
in 14]. 

In this paper we demonstrate that the disappearance 
of a neutron (or two neutrons) from its intra-nucleus 
state creates subsequent nuclear de-excitations that, for 
a few particular de-excitation modes, would lead to a 
unique experimental signature: a chain of time- and 
space- correlated events. Such a signature can be ob- 
served in the modern large-mass, low-background, and 
low-threshold detectors, e.g., in KamLAND 2]. Wc es- 
timate that the sensitivity of a search for such events 
in the KamLAND detector will result in improvement of 
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one-neutron and two-neutron disappearance limits by 3-4 
orders of magnitude beyond the existing lifetime limits. 

The idea that nucleon decay in the nucleus can be stud- 
ied by measuring particles and/or 7 rays accompanied by 
de-excitation of the nucleon hole was, to our knowledge, 
first considered by Y. Totsuka and H. Ejiri and 
used by the Kamiokande Collaboration Tt'I as discussed 
below. Also, recently the DAMA Collaboration 18] has 
searched in 6.5 kg of a liquid ^^^Xe scintillator detector 
operating in Gran Sasso National Laboratory for the de- 
cay products of the residual radioactive daughter nuclei 
left after the disappearance of one or two nucleons from 
a ^^^Xe nucleus and set new limits for two-nucleon dis- 
appearance. 

In sectional of this paper, we discuss the existing ex- 
perimental limits for n — > Sv, nn — *■ 2j/, and, in gen- 
eral, for disappearance of one or two neutrons. Our 
calculations for de-excitation branching ratios in intra- 
nucleus disappearance of one neutron from ^^C are dis- 
cussed in section IIIII for comparison, in section IIVI we 
also calculated the disappearance of a neutron from an 
^^O nucleus; section El deals with the calculations of de- 
excitation modes following the two-neutron disappear- 
ance from ^^C. The possible observation of some par- 
ticular de-excitation modes followed by the radioactive 
decay of daughter nuclei in the KamLAND detector, as 
well as the possible sources of background, are discussed 
in section IVII Our conclusions are presented in section 

IviTl 



II. EXISTING LIMITS FOR NEUTRON 
DISAPPEARANCE 

Experimental limits have been set for n — > decay 
in underground experiments [l9l |20| assuming that the 
whole Earth is the source of decaying neutrons. The 
underground detectors in this case measure the specific 
energy spectrum above the background level of atmo- 
spheric neutrino signals, resulting from decay-neutrinos 
interactions with the material of the detector. The best 
limit set from the analysis of the muon events induced 
by charge-current interactions in the detector |l9j| was 
T(n — > 3f^) > 5 • 10^^ years. The best limit from in- 
teractions was set in as T(n Sve) > 3 ■ 10^^ years. 
Due to the detection method, these limits are only valid 
for the specific decay of a neutron into particular types of 
neutrinos and cannot generally be applied to the neutron 
disappearance. 

A very inspirational approach has been used by the 
Kamiokande Collaboration who set the limit for 
n — > vi'D decay of r > 4.9 • 10^^ years, independent of the 
type of neutrinos emitted in the neutron decay. By virtue 
of the experimental method, this limit is also relevant 
for intra-nucleus neutron disappearance. The method is 
based on the detection of 7-rays, which could be residual 
de-excitation products of a highly excited ^^O nucleus 
left after the disappearance of a neutron from ^^O in the 



low-background energy range between 19 and 50 MeV. 
Because of the similarity to our approach, we are dis- 
cussing this method below in more detail. In this regard 
we also refer to the paper by H. Ejiri pjj. 

When a neutron disappears from the oxygen nucleus, 
the residual ^^O is often left in an excited state. In the 
nuclear shell model the excitation level is determined by 
the level that a neutron has inside the ^^O and by the 
width of this energy level (spreading width). In the ^^O 
nucleus, two neutrons occupy the lowest si/2 level, 4 neu- 
trons fill the next P3/2 shell, and the remaining 2 neutrons 
populate the highest pi/2 state. The disappearance of a 
neutron from the Pi/2 state would result in a minimal ex- 
citation where the residual ^^O with high probability will 
be left in the ground state and the difference of the bind- 
ing energy will be taken care of by undetected neutron 
decay products. The disappearance of a neutron from a 
lower state will create a hole that will result in a 
restructuring of the residual nucleus with possible emis- 
sion of 7 radiation. Thus, it was estimated in that 
de-excitation of the P3/2 neutron hole with a probability 
of 44% will produce 6.18 MeV 7-rays. The remaining 
part the ^^O nucleus will again be left in the ground state. 
When a disappearing neutron produces a hole in the Si/2 
shell, the excitation energy is large enough to exceed the 
separation energy for protons and neutrons in ^^O. In 
this case de-excitation proceeds mainly by emission of a 
proton, neutron, a, and, to a smaller extent, by emis- 
sion of 7s. Since the Kamiokande water-Cherenkov de- 
tector was not sensitive to low-energy protons, neutrons, 
as, and had a threshold for 7 and electron detection of 
7.5 MeV, the authors of paper concentrated only 
on the probability of de-excitation of the si/2-state hole 
via emission of 7s within an energy range of 19-50 MeV, 
where the experimental background was very small. This 
resulted in a rather small de-excitation branching ratio 
(Br) estimated in as (0.27-1.04) • 10""*. Due to this 
small branching ratio, the measured r/Br = 1.8 • 10'^^ 
years in the Kamiokande detector has resulted in t„^3^ 
disappearance limit of only > 5 • 10^^ years. 

For completeness, we need to mention that the Particle 
Data Group listings also include the n and 
n ^ lifetime limits of the J.-F. Glicenstein analysis 
[2l| r(n 3i/x) > 2.3-102'^ years andT(n bvx) > 1-7- 
10^^ years, based on the Kamiokande experimental data 
and on the idea that the disappearance of the neutron's 
magnetic moment should produce radiation. 

Existing lifetime experimental limits for nn —>■ 2v are 
even less stringent. The following two limits were deter- 
mined by the Frejus Collaboration j20|]: T{nn 2ve) > 
1.2-10^^ years and r(nn 2v^) > 6.0-10^^ years. These 
limits assume two-neutron decays into neutrinos and an- 
tineutrinos of e and /i flavors rather than the more gen- 
eral intra-nucleus disappearance of two neutrons. Re- 
cently the DAMA Collaboration looking for the de- 
excitation products of the residual radioactive daughter 
nuclei left after the two-neutron disappearance from a 
^^^Xe, has set the lower limit for nn — > 2^ (more gen- 
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erally for two-neutron disappearance) as r(rm —> 2v) > 
1.2 • 10^5 years. 



III. NEUTRON DISAPPEARANCE FROM '^C 

Large liquid scintillator detectors like KamLAND '5| 
and Borexino [l^l contain ^^C as a major component of 
the detector material. In the simple nuclear shell-model 
picture the ^^C ground state is described as a closed-shell 
nucleus i.e., for both protons and neutrons the lowest- 
lying Si/2 and shells are completely filled, while all 
higher (sub)shells are empty. This occupation of the nu- 
clear energy levels by neutrons is sketched in Figure ^ 
Thus, most of the neutrons available in liquid scintillator 
detectors for the baryon instability search are the neu- 
trons occupying Si/2 and energy levels in ^^C. 



A. Disappearance from the P3/2 shell 

If, in the simple shell-model picture, a neutron disap- 
pears from the shell in ^^C, the residual ^^C nucleus 
is always left in the ground state. However, in more so- 
phisticated shell-model calculations, it has been shown 
that the simple closed-shell picture of the ^^C ground 
state has to be remedied. Due to correlations, especially 
pairing effects, the pi/2 shell is also partially filled with 
protons and neutrons from the P3/2 level lying just 4.44 
MeV below (Fig^. Various shell- model calculations 
give partial occupation numbers of thepi/2 shell for both 
neutrons and protons in the range 0.75—0.8, correspond- 
ing to an average filling of the P3/2 shell between 3.2 
and 3.25. If we interpret the partial occupations by say- 
ing that approximately 60 percent of the time the P3/2 
shell is completely filled, while roughly 40 percent of the 
time two neutrons (the same for a proton pair) are in 
the Pi/2 shell, then there is a finite probability that a 
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FIG. 1: Occupation of energy levels by neutrons for the ^^C 
ground state in a simple shell model picture. One neutron is 
shown as disappeared from Si/2 level. 



neutron disappearing from the P3/2 state will leave the 
residual "'^^C nucleus in an excited state with spin and 
parity J'^ — 1/2^. There is only one state in -'^^C with 
= 1/2" below 10 MeV [Hi, which is the first ex- 
cited state at 2.0 MeV. Therefore, we assume that this 
is the only possible final state when a neutron dis- 
appears, while at the same time a pair of neutrons is in 
the pi/2 shell, i.e., we assume the corresponding spectro- 
scopic factor to be 1. Then the probability per P-shell 
neutron for this event is given by the product of 40 per- 
cent and 1/2, that is 20 percent, with an adopted error 
of ±5 %. The first excited state in ^-'^C will promptly 
de-excite by emission of a 7 with an energy of 2.0 MeV, 
and the ground state of "'^^C will subsequently undergo 
a j3'^ transition (the half-lifetime of ^^C is 20.4 min and 
Qec = 1-98 MeV). 

Thus, the disappearance of a neutron from the P3/2 
state will result, with a branching of 20±5% per each 
of four P-shell neutrons, in two time-correlated events: 
emission of monoenergetic 2.0 MeV 7 and (3^ emission 
with detectable energy in the range 1.02-1.98 MeV. 



B. Disappearance from the Si/2 shell 

A more complicated and rich picture emerges if a neu- 
tron disappears from the Si/2 state in ^^C. Excitation 
energy of the residual ^^C* in this case is very high; it 
exceeds the separation energy for proton, neutron, and a 
particles in ^^C, making them the primary emission prod- 
ucts of the highly excited nucleus. Qualitatively, emis- 
sion of a proton should be more frequent than emission 
of a neutron, since ^"'^C is a proton-rich nucleus. Pri- 
mary emission might not be sufficient for removing all 
the excitation energy and therefore might be followed by 
secondary de-excitation. Nuclear model calculations are 
needed in order to estimate these particular de-excitation 
branching ratios. 

The residual ^^C nucleus is left in a highly excited state 
with (assuming conservation of angular momentum and 
parity) total spin and parity J'" ^ ^ . The excitation 
energy E* of this state can be calculated by the dif- 
ference between the binding energy of the Si/2 neutron 
level in ^^C and the neutron separation energy of ^^C. 
While the latter is well known to be Sn = 18.72 MeV, 
no precise experimental or theoretical values exist for the 
binding energy of the Si/2 neutron level. In the case of 
the Sif2 proton level in ^^C the binding energy and the 
width have been measured in several electron scattering 
experiments ^^C(e, e'p) 253 and were determined to be 
Elg = 39 ± 1 MeV and F^'^, = 12 ± 3 MeV, respectively. 
Correcting E^g for the Coulomb shift of w 2.7 MeV for 
protons in -^^C gives us i?"g=41.7 ± 1 MeV for the bind- 
ing energy of the Si/2 neutron level in ^^C and, sub- 
tracting the neutron separation energy, we get the value 
i?* = 23± 1 MeV for the excitation energy in -'^^C. Deter- 
mination of width of the excited state caused by the dis- 
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appearance of Si/2 nucleon is not that straightforward. In 
terms of nuclear correlations that are responsible for the 
width of the excited state, the disappearance of nucleon a 
priori is not the same process as the knocking-out of a nu- 
cleon, since in the latter case a strongly-interacting par- 
ticle remains present in the final state inside the nuclear 
matter. Different approaches are possible here. Theoret- 
ical estimates for the spreading width of the si/2 hole in 
infinite nuclear matter [2^ gives a value of ~6 MeV. As 
another extreme one can use the spreading width of Sjuo 
state as determined in proton knock-out exeriments |25j | . 
We adopt what we think is a reasonable compromise here 
and assumed that the excited level at i?* = 23 ± 1 MeV 
in ^^C has a Lorentzian width of F = 7 MeV. At the end 
of the section we will show that our main results do not 
sensitively depend on the values adopted for E* and T in 
the whole range of variations mentioned above. 

In order to determine the de-excitation modes of the 
excited ^^C nucleus, i.e., the branching ratios in the var- 
ious particle emission channels and the energy spectra 
of the emitted particles, we used the statistical model 
of compound nuclear reactions (Hauser-Feshbach) [2^ . 
Within this model the branching ratio for the decay 
A* {c)Ac of an excited state of a nucleus A with energy E* 
and spin and parity to a final nucleus Ac by emission 
of a particle c is given by the ratio of the transmission 
coefficients: 



Brc{E*,r) 



Tc{E\r) 

Ttot{E*,J^) 



(1) 



Here Ttot is simply the sum of the transmission coeffi- 
cients of all decay channels open at E* : 



Ttot{E\J^) 



Tc{E\r) 



(2) 



which guarantees that the sum of the branching ratios 
is normalized to one. Each of the transmission coeffi- 
cients Tc is obtained by adding the contributions from 
all allowed final states in the residual nucleus: 

TciE*,.r) = J2TnE*,J-,E'!,Jlf,7:i!) + 



J ^ r,(£;*, j„7r,)p(£;„ (3) 



In this sum the decays to the ground state = 0) and all 
experimentally known excited levels {^i — 1,2,.. .^Nc) in 
the daughter- nucleus are explicitely taken into account. 
For excitation energies higher than the energy E^ of the 
last known level the sum is changed to an integration over 
the level density p up to the highest energetically allowed 
state at energy E* — Sc, where Sc is the channel sepa- 
ration energy for the ejected particle c. It can be seen, 
that the important ingredients of statistical model cal- 
culations are the particle and 7-transmission coefficients 



T and the level density of the excited states. Therefore, 
the quality of the calculated results within the statistical 
model is determined by the accuracy with which these 
components can be evaluated. 

The branching ratios and particle energy spectra we 
present in the following have been calculated with the 
statistical model code SMOKER f28j, which applies re- 
alistic optical potentials (reproducing nicely experimen- 
tal data) to determine the transmission coefficients and 
a sophisticated description for the level density (by a 
combination of a constant temperature formula and a 
back-shifted Fermi gas). Starting from an excited state 
with well-defined energy, angular momentum, and par- 
ity, the SMOKER program calculates the branching ra- 
tios for transitions via photon, neutron, proton, and a 
emission, which are the dominant de-excitation channels. 
Note that within SMOKER the transmission coefficients 
in Eq. (jS)) describing the transition from the excited state 
{E*, J^) to the state (E^, J^,tt^) in the daughter nucleus 
are correctly obtained by summing up the contributions 
from all quantum mechanically allowed partial waves: 



T^^iE*,r,Ei!,Jlf,7T^:) = 

J+S J^+Ja 

^ ^ Tc,SE*-Sc-En (4) 

Here angular momentum I and the channel spin s — 
jc + Jc couple to J = / -I- s. For the physical quantities 
typically needed in statistical model calculations (e.g., 
nuclear masses, level density parameters), SMOKER 
takes the experimental values, if they exist. Otherwise, 
theoretical models are applied. This is illustrated in Fig- 
ure m where the branching ratios for the decay of an 
excited state of ^^C* with spin and parity J'^ = 
into the photon, neutron, proton, and a channels are 
shown as functions of excitation energy. Note that all 
particle decay channels open at the correct thresholds 
{Sa = 7.54 MeV, Sp = 8.69 MeV, Sn = 13.12 MeV) 
because the experimental masses of the involved nuclei 
have been used. With increasing excitation energy, as 
soon as a hadronic decay channel opens, the branching 
for 7-decay (electro-magnetic interaction) becomes neg- 
ligible. While the contribution to the neutron channel 
rises immediately after the threshold, the branching into 
the proton channel needs 2-3 MeV above the threshold 
to become open because the Coulomb barrier has to be 
overcome. 

If particle emission decay leads to an excited level of 
the residual nucleus, we again calculate the branching ra- 
tios for the subsequent emission with the statistical model 
(iteration of SMOKER) . By keeping track of the energies 
of the ejected particles and photons during the cascade, 
and by weighting them with the corresponding branching 
ratios, the energy spectra for emitted neutrons, protons, 
as, and 7s can be determined. Since the amount of exci- 
tation energy in the system is significantly reduced with 
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FIG. 2: Energy dependence of de-excitation branching of a J'^ = ^ ^^C* state by emission of a 7, neutron, proton, or a 
particle as built into the statistical model (SMOKER code). 



every particle emission and eventually drops below the 
particle emission thresholds, we found that the cascade 
could practically be terminated after two emission steps, 
i.e., the residual daughter nuclei after two iterations were 
assumed to be in their ground state. 

The SMOKER code has been applied extensively to 
many astrophysical problems, mainly where nuclei with 
mass number ^ > 20 are involved. In such a case the 
requirement of a high level density, which is one of the 
main assumptions of the statistical model, is usually met, 
and the calculated cross sections are generally in good 
agreement with the experimental data [23 ■ For the light 
nuclei considered here, an application of the statistical 
method might seem problematic because the level den- 
sity in light nuclei can be small, and in many reactions 
with light nuclei, significant direct contributions have 
been found. However, by studying the de-excitations of 
^^C and ^^O nuclei that had been excited by inelastic 
neutrino scattering, it was found that the branching ra- 
tios for decay into the proton and neutron channels, cal- 
culated within the statistical model, were in very good 
agreement with the branching ratios obtained from a con- 
tinuum random phase approximation |30j |. That means 
that at certain conditions, direct contributions can be ef- 
fectively accounted for in the statistical model. Although 
at specific excitation energies, the branchings calculated 
within both models deviated by a factor of 2, the proton 
and neutron branchings averaged over various excitation 
spectra corresponding to different neutrino energies were 
found to match within 1-17%. As for our calculations 



in this paper, due to the spreading width of the excited 
state, the branching ratios from the statistical model are 
also averaged over a range of excitation energies. There- 
fore, we conservatively adopt the relative error in the 
branchings for one- and two-step decays to be 20% and 
30%, respectively. 

We now turn to the results of our calculations of the 
final states of nuclear de-excitations following the disap- 
pearance of a neutron from the si/2 shell in ^^C. Results 
of these calculations are shown in Table |l| They were 
obtained for an excited J'" = level at E* = 23 MeV 
in ^^C with an associated Lorentzian width of r=7 MeV. 

The first upper block in the table just shows the 
branching ratios for one-step emission into the 7, neu- 
tron, proton, and a channels. In the second block the 
branching of the first-step neutron channel is divided 
into the contributions directly to the ground state and 
the contributions to the excited states in the daughter 
nucleus ^"C, which undergo subsequent second-step 7, 
neutron, proton, or a emission. In the same way, the 
first-step proton and a channels are split up in the lower 
blocks of Table In comparison to the time scale rel- 
evant for the detection of the ejectiles, the emission of 
de-excitation particles occurs simultaneously. Hence the 
"symmetric" channels like "'^"'^C(n,p) and ^^C(p,n) will 
lead to identical signatures in the detector and therefore 
we list the sum of these branchings in the lowest block of 
TableQ] Next to the two-step de-excitation mode, shown 
in the first column of Table the residual daughter nu- 
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TABLE I: Branching ratios for ^^C* de-excitation modes 
after a neutron disappears from the Si/2-state in '^^C and 
experimental signature (number of time-correlated hits) for 
observation of these modes in a large liquid scintillator detec- 
tor. For the final states marked with * , due to the very large 
lifetime of the daughter nuclei, the experimental signature is 
assumed to be a single hit. 



Decay 


Daughter 


Mode 


Exp. 


mode 


(decay, T1/2 or T, Qec or Q^-) 


% 


sign. 








(hits) 




''Cgs(/3+; 20.4 m, 1.98 MeV) 


0.7 


2 


"C(n...) 


i 


13.8 


i 


"C(p...) 


i 


64.4 


i 


"C(a...) 


i 


21.1 


I 


''C(n) 


"'Cgs(/3+; 19.3 s, 3.65 MeV) 


3.0 


3 


"C(n,7) 


'''Cgs(/?+; 19.3 s, 3.65 MeV) 


2.8 


3 


"C(n,n) 


^C{I3+; 0.127 s, 16.5 MeV) 


0.06 


4 


"C(n,p) 


''B(p+2a; 0.54 keV, 1.07 MeV) 


5.7 


2 


"C(n,a) 


'*Be(2p+a; 92 keV, 4.3 MeV) 


2.2 


2 


''C(p) 


'"B(stable) 


2.9 


1 


"C(p,7) 


^''B(stable) 


19.0 


1 


"C(p,n) 


''B(p-f2Q; 0.54 keV, 1.07 MeV) 


1.4 


2 


"C(p,p) 


^Be(stable) 


7.2 


1 


"C(p,a) 


<^Li(stable) 


33.9 


1 


-C(a) 


'Be{l3+; 53.3d, 0.86 MeV) 


5.2 


1* 


"C(a,7) 


'■Be(/3+; 53.3d, 0.86 MeV) 


4.2 


r 


"C(Q,n) 


^Be(2p+a; 92 keV, 4.3 MeV) 


0.3 


2 


"C(a,p) 


8Li(stable) 


3.4 


1 


"C(Q,a) 


^H(/3-; 12.3y, 18.6 keV) 


8.0 


1* 


^'C{n,p} 


''B(p+2q; 0.54 keV, 1.07 MeV) 


7.1 


2 


"C{n,a} 


^Be(2p+a; 92 keV, 4.3 MeV) 


2.5 


2 


"C{p,a} 


'^Li(stablc) 


37.3 


1 



clei are listed in the second column. Some of these nuclei 
are radioactive and can provide a time-correlated signal 
with the emitted particles in the detection process. The 
third column gives branching ratios in percentages per 
one neutron in the si/2 state. 

Energy spectra for de-excitation 7s and neutrons are 
shown in Figures |21 and 0] respectively. Figure |21 shows 
the calculated 7-spectra for all de-excitation modes (dot- 
ted Hne), for ^^C*{'yy^Cgs transitions (dashed line), and 
for the transition ^^C* {n,^)^°Ggs (solid line). The last 
spectrum is dominated by a strong monoenergetic line 
at 3.35 MeV, corresponding to the decay of the first ex- 
cited 2+ state in ^°C. Figure 01 shows the energy spec- 
tra for the emitted neutrons. While the solid line de- 
picts the summed neutron spectra from all processes fol- 
lowing the disappearance of a Si/2 neutron, the dotted 
line just shows the neutron energy distribution from the 
^^C*(n, 7)-'^°Cgs mode, and the dashed line relates to the 
neutrons from the decay ^^C*(n)^°Cgs. The spectrum 
of neutrons from the "'^^C*(n, 7)^''Cgs mode reflects the 
width of the initial neutron level in ^^C. 

In liquid scintillator detectors the emission of de- 
excitation 7s, protons, and as will be seen as a prompt 



TABLE IL Branching ratio (in percent) into the 
^^C* (n,7)^'^Cgs decay channel as a function of binding energy 
and spreading width of Si/2 neutron level in ^^C. 
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signal, and we assumed here that these particles can- 
not be distinguished. The neutrons, although emitted 
at the same time as other de-excitation products, will 
result in a double signal. First, they will quickly slow 
down by elastic collisions with hydrogen (which, besides 
^^C, is another essential component of the liquid scintil- 
lator) and thus will also contribute to the prompt com- 
ponent of the signal. Second, thermalized neutrons will 
diffuse in the liquid scintillator for typically 200 /isec be- 
fore they are captured by hydrogen atoms with emis- 
sion of a detectable 2.2 MeV photon. Thus, neutrons 
can provide two signals in the chain of detected events. 
Residual radioactive nuclei left after one- or two-particle 
de-excitations will provide another signal in the chain of 
events initiated by the disappearance of the neutron, all 
coming from the same point in the detector within the 
reconstruction space resolution and within certain corre- 
lation time. We also note that positrons from f3~^ decay 
of the daughter nuclei in Tabled can be seen in the liquid 
scintillator detector with efficiency close to 100%. 

The chain of time- and space-correlated events creates 
a unique signature for the detection of neutron disappear- 
ance in low-threshold liquid scintillator detectors. The 
right-most column of Tabled shows the experimental sig- 
natures of the corresponding decay modes in terms of 
time- and space-correlated hits in the detector. One can 
see that the decay chains: ^^C* n + j+^^Cgs, with 
a branching of 2.8%, ^^C* — > n+^°Ggs, with a branch- 
ing of 3.0%, and ^-^C* — > n -f n+^Cgs, with a branch- 
ing of 0.06% provide 3-hit or 4-hit signatures. Therefore 
these decay modes should be very immune to possible 
accidental-coincidence background. 

As the initial values for the excitation energy E* and 
spreading width F are not precisely known, in order to 
estimate the corresponding systematic stability of our re- 
sults, we varied these quantities within their uncertainty 
range: ±1 MeV for E* and from 6 MeV to 15 MeV for 
F. Exemplary for all decay modes, the effect of the varia- 
tions on the results are shown in Table ITll for the branch- 
ing to the ^^C* ^ n + J -|-^° Cgs channel, which we con- 
sider to deliver the most robust signature for a neutron 
disappearance. 

One can see from Table Ull that the predicted branch- 
ing ratio is rather insensitive to the variation of E* and 
F. Since, as mentioned above, an error of 30% should be 
assigned to the branching ratios, we consider the system- 
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FIG. 3: Spectra of 7s from de-excitations of si/2 hole of neutron disappearance in ^^C. The dotted line is for all 7s from all 
de-excitation processes; the dashed line is for ^^C* (7)^^053 transition, and the solid line is for the de-excitation mode n + 7 
with a residual ^°Ggs nucleus. The energy bin width is 0.2 MeV. 
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FIG. 4: Spectra of neutrons from do-oxcitations of Si/2 hole of neutron disappearance in ^'^C. The solid line is for all neutrons 
from all de-excitation processes; the dotted line is for the process with emission of de-excitation n, 7, and ^'^Cgs residual nucleus; 
the dashed line is for the process with de-excitation via emission of neutron and ^"Cgs residual nucleus. The energy bin width 
is 0.2 MeV. 
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atic error coming from the uncertainty in E* and T to be 
negligible. 

In this paper we assume that decay products on neu- 
tron disappearance (e.g. n wO ) have total momen- 
tum equal to the momentum of the initial neutron. That 
can be justified by the assumption that disappearance 
of a neutron occurs at very short distances and for very 
short times relative to the distances and times that are 
characteristic of nuclear processes. Therefore, one can 
expect that the momentum of the disappearing neutron 
is taken away by the decay products except possibly for 
a small recoil given to the residual nucleus as a whole. 
Our main result, the branching ratios for the residual 
nucleus de-excitation, will not be affected by the recoil 
motion. However, spectra of emitted protons, as, and 
neutrons might be smeared out by this motion. From 
experimental point of view this smearing will not lead 
to any principal effects, it might only slightly increase the 
uncertainty of the efficiency of detection of these particles 
due to the finite visible energy threshold. 

IV. NEUTRON DISAPPEARANCE FROM ^'^O 

Neutron disappearance from the Sx/2 state in ^^O fol- 
lowed by 7 de-excitation from a highly excited state has 
been considered in |17| for a Kamiokande experiment 
with a water-Cerenkov detector. The branching ratio for 
this de-excitation mode was calculated in this paper. It 
is interesting to compare our branching ratio calculations 
for the initial ^^C nucleus with this result. Integrating 7 
spectrum of ^^C* de-excitation, which is shown in Figure 
13 as a dotted line, within the energy range from 19 to 
50 MeV, as suggested by Ref. found a branch- 

ing ratio of 1.1-10~*. This is similar to the si/2 hole 7 
de-excitation branching ratio in ^^O, calculated in ^l_7j as 
(0.27 - 1.04)-10-1 

We have also calculated the spectrum of final de- 
excitation states for an excited ^^O nucleus, originated 
by a neutron hole in the S1/2 neutron level in ^^O. These 
branching ratios are shown in Table UTTl They were ob- 
tained by SMOKER code calculations assuming the val- 
ues E* = 29 MeV and F = 7 MeV respectively for the 
energy and spreading width of the excited state in ^^O. 

Figure [3 shows in a logarithmic scale the spectrum of 
7s from Si/2 de-excitations of ^^O*. For the 7 branching 
ratio, integrated over the energy range 19—50 MeV, we 
obtain the value 1.4-10"''. By varying the energy level 
and width of the 5^/2 state in ^^O by ± 1 MeV, we found 
the branching to change in the range (1.3—1.5) •10"'*. 
Furthermore, taking into account the uncertainty in the 
threshold energy of 19 MeV due to the energy scale cali- 
bration of ± 0.6 MeV, the range becomes (1.2— 1.6)-10"'*. 
Finally, estimating the systematic uncertainties of the 
SMOKER code, we conclude that 7 branching for de- 
excitation of the Si/2 hole in -'^O within the energy range 
19-50 MeV is (1.4 ± 0.7) • 10"^, which is in good agree- 
ment with the result (0.27— 1.04)-10~^ of the calculations 



TABLE III: Branching ratios for ^^O* de-excitation after neu- 
tron disappearance from a s-1/2 hole state in ^®0. 
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from Ref. [T3|- 



V. TWO-NEUTRON DISAPPEARANCE 
FROM ^^C 

Since the disappearance of two neutrons would be a 
process occurring at a very short distance, we should 
consider both disappearing neutrons to be in the same 
subshell, and thus we will neglect the cases when one S- 
shell and one P-shcU neutrons are disappearing together. 



A. Disappearance from the pi/2 or P3/2 shell 

In the simple shell-model picture, a disappearance of 
two neutrons from the P3/2 shell would always lead to 
= 0+ ground state of ^"0. However, due to corre- 
lations and particularly due to nucleon pairing, as was 
discussed earlier, there is a probability of ~ 40 % of find- 
ing a pair of neutrons in the pi ^2 level. We again assume 
that the disappearance of two neutrons would occur due 
to a very short-range interaction, and therefore expect 
the disappearance of paired neutrons in either a Pi/2 or 
a P3/2 state to be more likely than the disappearance of 
a pair of one Pi/2 and one P3/2 neutrons. Nevertheless, 
we will discuss both cases below. 
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FIG. 5: Spectrum of 7s from de-excitation of a Si/2 hole resulting from neutron disappearance in O. The energy bin width 
is 0.2 MeV. 



The disappearance of paired neutrons with total spin 
would imply certain selection rules for the final state of 
disappearance. Thus, from the conservation of angular 
momentum, the transition nn — > vD with A{B — L) = —2 
is possible while transitions rin — > w or nn vD are 
not. The disappearance of two neutrons from the pi/2 
level will lead to the ground state of ^"C. The disappear- 
ance of two neutrons from a state, while two other 
neutrons are in a Pi/2 state, might result in the excited 
state of "'^^C with J'^ = 0"''. Since there are no excited 
states of ^"C with this spin and parity below the pro- 
ton separation threshold, the residual de-excitation will 
be dominated by i^C* p-K^B, and will quickly de- 
cay into p + 2a. Such de-excitation will lead only to a 
single-hit signature in the liquid scintillator detector and 
therefore is not favorable for detection. 

A disappearance of one P3/2 and one pi/2 neutron 
would lead to an excited state in ^"^C with spin and parity 
J'^ — 2+ or 1+. There is not any low lying energy level 
known in ^"C with spin and parity 1^ 24j. However, 
there are 3 low lying energy levels stated in -'^^C (two 
are well established and one is assumed) with spin and 
parity 2+ that could be populated after 2-neutron 
disappearance. To our knowledge, no experimental or 
theoretical information on the spectroscopic factors for 
these possible final states exists here. However, as 2 out 
of the three excited 2+ states in lie above the proton 
or even a separation threshold, their decay will proceed 
by emission of a low energetic proton or a-particle, lead- 
ing to the ^B and ^Be ground states, respectively. As ^B 



immediately decays into a proton and two a particles, 
and the ^Be into 2 protons and a, this will lead only to 
a single-hit signature in the liquid scintillator detector. 
Only the population of the first excited 2+ state at 3.35 
MeV will produce a space and time correlated two-hit 
event in the detector. This state decays by emission of a 
photon with energy 3.35 MeV, followed by the /S"*" decay 
of the residual ^°Cgs with a correlation time of 27.8 s and 
detectable energy between 1.74 and 3.65 MeV. As an up- 
per limit for the probability of this event we adopt the 
value 40 percent x 4/6 (probability to pick a (p3/2,Pi/2) 
pair from two P3/2 and two pi/2 neutrons) / 3 (number 
of low lying 2+ states in ^'^C) or 9 % per ^^C atom, with 
an estimated error of ± 5 %. 



B. Disappearance from the Si/2 shell 

For the disappearance of two neutrons from the Si/2 
state, we estimated the excitation energy of i^C* in the 
following way: bound energy of two neutrons in the Si/2 
state is assumed as 78 MeV; separation energy of two 
neutrons from ^^C is 31.8 MeV; thus the excitation en- 
ergy is ^46 MeV. With respect to the spreading width 
of a double Si/2 hole state one might be inclined to as- 
sume that it should be larger than the width of a sin- 
gle Si/2 hole (7 MeV). However, essential contribution 
to the spreading width of the single Si/2 hole is due to 
the pairing correlations that should be significantly re- 
duced for a double Si/2 hole state. Therefore, a width of 
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r=7 MeV has been assumed for a double si/2 hole state. 
The results of SMOKER code calculations for this case 
are shown in Table IIVI and were found to change their 
relative value by less than 3.5 percent when varying the 
spreading width between r=4 and 15 MeV. 

There are two final-state modes in this table that would 
result in a 3-hit event signature in the liquid scintillator 
detector, each with a branching of about 6% per S-shell 
neutron. Figure El shows the energy spectra of neutrons 
for these de-excitation modes. Note that both spectra 
reflect an initial Lorentz-shape distribution of excitation 
energy. For the mode ^°C* n + p+^Bgs, the proton, 
which comes from the decay of the first excited state in 
^C, is monoenergetic with energy 0.922 MeV. 



VI. DISCUSSION 

The KamLAND detector P| is perfectly suited to 
search for the neutron disappearance processes in car- 
bon, i.e., for the observation of sequences of nuclear de- 
excitation events. The active mass of the KamLAND 
detector filled with a liquid scintillator (formula ^ CH2) 
is ^ 900 tons; the trigger threshold can be as low as ~ 0.2 
MeV; the photon energy resolution at 1 MeV is < 8%. 
Good time and energy resolution allow reconstruction of 



TABLE IV: Branching ratios for C* de-excitation after two- 
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neutron disappearance from Si/2-state of ^^C. 
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an event position in the detector with the accuracy down 
to several cm. Low background in the KamLAND experi- 
ment is due to a deep underground location (2,700 meters 
of water equivalent), water shield, buffer oil shield, veto 
system, radon shields, and liquid scintillator purification 
that provide concentration of ^^^U, ^^^Th, and '*'^K at 
the level below IQ-^^ g/g |29l |. 



Experimental signatures 

We have calculated in sections IIIII and above the 
branching ratios for various de-excitations following the 
disappearance of one or two neutrons from the ^^C nu- 
cleus. Not all de-excitation modes are favorable for de- 
tection in the experiment, but some of them form the 
chains of time- and space-correlated events in the detec- 
tor. Thus, the process ^^C* n + j+^'^Cgs from the 
disappearance of a neutron from the si/2 state will look 
like a 3-hit event in the detector. The first hit in the 
detector will be produced by monoenergetic 7s with an 
energy of 3.35 MeV corresponding to the de-excitation of 
a 2+ level of and by slow-down interactions of neu- 
tron in the scintillator. The first hit will be followed by a 
second hit, which is due to neutron capture by hydrogen 
in the scintillator (correlation time is ^ 200/xsec). The 
third hit will be due to delayed /J"*" decay of ^"C^s with 
correlation time 27.8 sec and detectable energy in the 
range 1.74-3.65 MeV. The first two hits will look similar 
to a signature of the reactor antineutrino interactions: 
positron plus delayed captured neutron. The detection 
of reactor antineutrinos is the main purpose of the Kam- 
LAND experiment . The expected number of reactor 
antineutrino interactions in KamLAND in the absence 
of neutrino oscillations is ^ 600 per kt per year with the 
background less than few percent . Using the measure- 
ments performed in the KamLAND Collaboration [2^ . 
one can conservatively estimate the background for an 
accidental 3-hit coincidence (first hit with a threshold of 
2.8 MeV, second hit with a threshold of 1.6 MeV; win- 
dow for first two hits is 0.54 msec; threshold for the third 
hit is 2 MeV; time window for the third hit is 60 sec; 
position-correlation cut of Im) to be below one event 
per year. This accidental background rate can be mea- 
sured in the detector and subtracted. Using branching 
for ^-'^C* n + 'j+^'^Cgs from Tableland assuming no 
other sources of events with such a signature (e.g. from 
atmospheric neutrino interactions), for one year of ex- 
position time in the KamLAND detector, one can, with 
90% CL, set a limit for n vvv decay at the level 5T0^^ 
years. This is 1,000 times better than the present limit 
for this decay mode fl|. 

Detection of the 3-hit process ^^C* — * n+^'^Cgs will be 
similar to the one discussed above, probably with a differ- 
ent threshold for the first hit that will be formed by the 
slow-down interactions of neutrons in the liquid scintil- 



lator. Detection of the 4-hit process C 



would be a spectacular observation; unfortunately its 
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FIG. 6: Energy spectrum of neutrons from de-excitations of a Si/2 hole resulting from two-neutron disappearance in ^^C. The 
solid line is for transition ^°C* — > n+^Cgs and the dotted line is for ^''C* — > n + p+^'Bgs. The energy bin width is 0.2 MeV. 



branching, as follows from Tabled is much lower than for 
two processes with the 3-hit signature discussed above. 



Non-accidental background 



The two- hit events listed in Table H] should be divided 
into two categories. Events with one of the de-excitation 
particles being a neutron belong to the first category; in 
KamLAND these events will have a signature similar to 
that for anti-neutrino reactor events with a correspond- 
ing background and will be difficult to disentangle. The 
second category should include events with radioactive 
daughter nuclei. Events resulting from the neutron dis- 
appearance from the P3/2 state of ^^C should also belong 
in this category. Accidental coincidence background will 
be more of a problem for the events with long lifetime 
daughter isotopes therefore only modes with a production 
of ^^Cgs might have practical interest for the detection. 
Correlation time between two hits in this case is 20.4 min 
due to lifetime of ^^C^s. The statistical fluctuations of 
the background will determine the sensitivity limit for a 
corresponding process search. Optimistically, this mea- 
surement alone can improve the existing lifetime limit for 
n vvD by approximately an order of magnitude. 



De-excitation final states originated by the disappear- 
ance of two neutrons from ^^C provide two interesting 3- 



hit signatures: ^^^C* 
and i°C* - 



Cgswith a branching of 6.2% 



n+p+°B with a branching of 6.0%. The rel- 
atively short lifetime of the daughter nuclide makes such 
events free of accidental background. 



It is interesting to see whether atmospheric neutrino 
events can provide a signature similar to ^^C* n + 
7-|-^°Cgs, ^^C* n-|-^°Cgs, or other processes. We an- 
ticipate in the future to study this interesting question in 
a separate paper treating nuclear de-excitations produced 
by atmospheric neutrinos and anti-neutrinos in "'^^C. Such 
a study will also be interesting from the point of view of 
identification of atmospheric neutrino and anti-neutrino 
interactions in liquid scintillator detectors through the 
detection of associated nuclear chains of de-excitations 
and radio-nuclide decays. In the respect of the back- 
ground for one-neutron or two-neutron disappearance ad- 
dressed in this paper, few qualitative comments can be 
made here. Charged current atmospheric neutrino reac- 
tions have a rate of about 100 interactions per kt-detector 
per year. Due to the relatively high energy spectrum 
of atmospheric neutrinos, leptons produced in charged 
current interactions usually have an energy larger than 
the energy typical for nucleus de-excitation. We remind 
the reader that in the Kamiokande experiment T.7j the 
energy range between 19—50 MeV was practically free 
of background. Thus, one can hope that charge cur- 
rent processes initiated by atmospheric neutrinos will 
not produce significant background for neutron disap- 
pearance processes, i.e., the chain of events with the first 
signal in a few-MeV range. However, detailed simulations 
of atmospheric neutrino interactions with ^^C using the 
SMOKER code for treatment of nuclear de-excitations 
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are needed to understand this background quantitatively. 
In such simulations, the production of low energy charged 
leptons has to be estimated together with de-excitation 
of ^^N* and ^^B* nuclei with a continuous spectrum of 
excitations. Neutral current atmospheric neutrino inter- 
actions should be also followed by de-excitations of ^^C* 
(see paper (sij) and might, for example, potentially lead 
to the process ^^C* ^ n + n + j+^'^Ggs with one of the 
neutrons lost at the edge of the fiducial volume of the 
detector. 



VII. CONCLUSIONS 

In this paper we show that modern high-mass, low- 
threshold, and low-background liquid scintillator detec- 
tors can observe unique signatures of chains of time- and 
space-correlated events, formed by nuclear de-excitations 
following the neutron or two-neutron disappearance from 
the nucleus. The KamLAND detector 2j is perfectly 
suited to search for such processes originating from ^^C 
nuclei. The sensitivity of search for a neutron disap- 
pearance (e.g. n vvv ) or two-neutron disappearance 
(e.g. nn — > vv) in the KamLAND detector can be ex- 
tended by 3-4 orders of magnitude relative to the present 
lifetime limits for these processes. Non-accidental back- 



ground for neutron and two-neutron disappearance from 
atmospheric neutrino interactions, although expected to 
be small, needs to be calculated in future work. Assum- 
ing that nucleons can decay only into known particles 
with conservation of electric charge, 4-momentum, and 
angular momentum, the search for intra-nucleus neutron 
and two-neutron disappearance in KamLAND can lead 
either to the discovery of long-awaited baryon instability 
or to the improvement of the mode-independent nucleon 
decay lifetime limit from the present 1.6-10^^ years to or 
above ~10'^° years. 

We think that the calculations performed in this paper 
with regard to the KamLAND detector and ""^^C nuclei 
can be extended to other high-mass, low-threshold, and 
low-background detectors using different nuclei. We plan 
to continue exploring these possibilities. 
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